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ABSTRACT 

Peculiar velocities are an important probe of the growth rate of mass density fluctuations in the 
Universe. Most previous studies have focussed exclusively on measuring peculiar velocities at inter- 
mediate (0.2 < z < 1) redshifts using statistical redshift-space distortions. Here we emphasize the 
power of peculiar velocities obtained directly from distance measurements at low redshift (z < 0.05), 
and show that these data break the usual degeneracies in the fl my o~ o's.o parameter space. Using only 
peculiar velocity data, we find S7 m o = 0.259 ± 0.045 and as,o = 0.748 ± 0.035. Fixing the amplitude 
of fluctuations at very high redshift using observations of the Cosmic Microwave Background (CMB), 
the same data can be used to constrain the growth index 7, with the strongest constraints coming 
from peculiar velocity measurements in the nearby Universe. We find 7 = 0.619 ± 0.054, consistent 
with ACDM. Current peculiar velocity data already strongly constrain modified gravity models, and 
will be a powerful test as data accumulate. 



o 
u 

Oh 
I 

o 

> 

00 

en 

o 



X 

S-H 



1. INTRODUCTION 

In the standard cosmological model, the Universe is 
dominated by cold dark matter combined with a cos- 
mological constant or dark energy (ACDM). While the 
existence of dark matter is supported by dynamical tests, 
gravitational lensing and fluctuations in the Cosmic Mi- 
crowave Background (CMB), the evidence for a cos- 
mological constant is primarily through geometric con- 
straints on the redshift-distance relation (standard can- 
dles, standard rulers) and he nce on the expansion h istory 
of the Universe (for example Komatsu et al.|[2011| here- 
after WMAP7-|-BAO-|--ffo)- Modified gravity theories, 
howeve r, can mimic t he expansion history of the ACDM 
model. Linder ( 2005 ) has emphasized that it is essential 
to measure the growth of structure as a function of cos- 
mic time as this allows one to break this degeneracy. He 
also shows that for many models, the logarithmic deriva- 
tive of the growth of structure can be parametrized as 



a hi a 



(1) 



where z is the redshift, g is the linear perturbation 
growth factor, a = 1 /(1 + z) is the expansion fa ctor and 



7 is 0.55 for ACDM ( |Wang fc Steinh ardt 199 8J. In con 
trast, for example, 7 = 0.68 in the |Dvah" et al.| ( 2000 



hereafter PGP) br aneworld modified gravity model! Lin- 
der fe Cahn|[2007| 



There are several ways to measure the amplitude of 
the dark matter power spectrum at redshifts lower than 
that of the CMB, including cosmic shear from weak 
gravitational lensing and the abundance of rich clus- 
ters. Another promising way to pro be the growth rate 
of structure is via peculiar veloc ities ( Guzzo et al.||2008~ 



There are two ways to measure peculiar velocities. The 
first method is statistical: given a galaxy redshift survey, 
the distortion of the power spectrum or correlation func- 
tion in redshift space depends on (3 = f/b, where b is a 
galaxy bias parameter ( Kaiser|1987 ). On large scales, we 
may assume that linear biasing holds, i.e. b = a St g/as, 
where as is the root-mean-square density contrast within 
an 8 Mpc/ft, sphere, h is the Hubble parameter in units of 
100 km/s/Mpc, and the subscript "g" indicates the fluc- 
tuations in the galaxy density, whereas no subscript in- 
dicates fluctuations in the mass density contrast. Galaxy 
redshift surveys also allow one to measure as, g directly, so 
one can combine the observables to obtain the combina- 
tion fas = (3&8,g- By combining redshift-space distortion 
(RSD) measurements of fas at different redshifts, one 
can study the growth of linear structures over a range of 



redshifts (Guz zo et aLJ20 08; Son gfc Percival|2009| [Blake 
|et al.||2011b[ |Samushia et al'||2oft| ^ 

A second method is to measure peculiar velocities di- 
rectly, by measuring distances to individual galaxies (via 
standard candles or standard rulers), and comparing 
these distances to their redshifts. We refer to this method 
as "measured distance" (MD). The measured peculiar 
velocities can then be compared with the peculiar veloc- 
ities predicted from a galaxy density field, S g (r), derived 
from an independent redshift survey, under the assump- 
tion that fluctuations in the mass are linearly re lated to 
those in the galaxy density. Specifically, we have ( Peebles 



1993) 



v(r) 
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Kosowsky & Bhattacharya 2009). Pec uliar velocities are This comparison of the two measured quantities v(r) and 
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directly proportional to the derivative of the growth fac 
tor, i.e. proportional to / 
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S g (r) yields a measurement of the degenerate combina- 
tion = f/b, which can be converted into the combina- 
tion fas in the same way as for RSDs. Note that because 
this is a velocity-density cross-correlation, unlike RSDs, 
it is not affected by cosmic variance. Instead, the noise 
arises from the precision of the MDs themselves. 
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Table 1 






Measurements 


of /<rg from the literature, 


and 






derived 7 






Label 


z 


f&8 


7 


Ref 


THF 


0.02 


0.398 ± 0.065 


56+0' 11 
U - OD -0.09 


1 


DNM 


0.02 


0.314 ±0.048 


71+0 .10 
u -' 1 -0.09 


2 


6dF 


0.07 


0.423 ±0.055 


54+0-09 
UiO ^-0.08 


3 


2dF 


0.17 


0.510 ±0.060 


43+0- 10 
u -*°-0.08 


4 


LRG1 


0.25 


0.351 ± 0.058 


77+0- 16 
u -' '-0.14 


5 


LRG2 


0.37 


0.460 ± 0.038 


55+009 


5 


WZ1 


0.22 


0.390 ± 0.078 


67+0 .19 
u -°'-0.15 


6 


WZ2 


0.41 


0.428 ± 0.044 


64+0 .12 
U - Dqi -0.11 


6 


WZ3 


0.6 


0.403 ± 0.036 


76+0 .14 
U -' D -0.13 


6 


WZ4 


0.78 


0.493 ± 0.065 


0.38«|* 


6 


BOSS 


0.57 


0.415 ±0.034 


71+°' 12 
u -' 1 -o.n 


7 


VVDS 


0.77 


0.490 ±0.180 


40+0' 89 
u '* u -0.59 


8 



References - (1) |Turnbull et al.l (|2012 

(2) Davis et al. (2011); (3) Bcutlcr et al. (2012 
(4) Song & Fcrcival (2009); (5) Samushia et a] 



(2012); (6) Blake et al. (2011a[| ; (7) I Reid et al 
l 2DT2 i; (8)|Ouzzo et al. pU08); Song & Fercival 
(2009) 



Note. - - The ACDM distance-redshift rela- 
tions are assumed when calculating the AP dis- 
tortion, but the growth factor /eg is free. Derived 
7 also uses WMAP measurement of fluctuations 
at z CMB 

The two peculiar velocity probes are complementary: 
RSDs require large volumes, driving one to surveys at 
higher redshifts. MDs have errors which are a constant 
fraction of distance. Hence the error in peculiar velocity 
in units of km s _1 increases linearly with distance and 
so MD surveys are necessarily restricted to low redshifts. 
However, as we will show it is the lowest redshift data 
that have the most "lever arm" for constraining the cos- 
mological parameters considered here. The important 
point is that by combining high and low redshift mea- 
surements of f(z)ag(z), we can break many degeneracies 
in the cosmological parameters. 

An outline of this Letter is as follows. In Section 2, we 
present the data used in our analysis. Section 3 presents 
fits in the r^o-oxo plane for the ACDM class of models. 
In Section 4, we allow 7 to vary, but use CMB data to 
constrain the amplitude of cts^cmb)- We discuss future 
prospects in Section 5 and summarize in Section 6. 

2. DATA 

We will combine measurements of fa s from the two 
distinct methods discussed above. The majority of the 
data are from RSD me asurements , many o f which were 
previously compiled in Blak e et al. (1 2 01 lb ), but adding 
recent results from B OSS ( |Reid et al.||2012| ) and 6dFGS 
( Beutler et al.|[2012 ). These are summarized in Table 1. 

As discussed above, it is also possible to measure 
fas from the comparison of peculiar velocity and den- 
sity fields, where peculiar velocities are obtained directly 
from MDs. Most of the results from thi s method (up 
to 2005) w ere summarized in Table 3 of Pike & Hud 



son 



( 2005 hereafter PH) . They showed that results from 



estimates of fag were consistent with 0.41 ±0.03. How- 
ever, since both the peculiar velocity and redshift survey 

2 After correcting for their assumed 7 = 0.6 and making an 



datasets overlap somewhat, the measurements are not 
independent and it is therefore difficult to assess the un- 
certainties. 

Its possible to obtain a value of fas using non- 
overlapping datasets: by comparing peculiar velocit ies 
from the "Composite" sample of |Watkins et al. ( 2009 ) to 
predictions derived from the IRAS Point Sourc e Catalog 
Reds hift (P SCz) survey galaxy density field ( |Saunders| 
et al.||2000[ ), we find A = 0.49 ± 0.04 after marginalizing 
over residual bulk flows. This yields fas — 0.37±0.04, in 
good agreement with the PH mean value quoted above. 
These uncertainties may underestimate possible system- 
atics, however, since only one density field (PSCz) is used 
and, furthermore, the same comparison method is used 
for all peculiar velocity samples. 

Therefore, in order to make a fair (albeit conservative) 
estimate of the uncertainties in the methods, in this Let- 
ter we focus on only two recent measurements of fag 
from MP su rveys. The first of these is from |TurnbulI| 
et al.| ( |2012[ hereafter THF), who compiled the 'First 
Amendment' (hereafter Al) set of 245 peculiar veloc- 
ities from Type la supernovae (SNe) with z < 0.067. 
The Al peculiar velocities were compared to predic- 
tions from the IRAS PSCz galaxy density field, yielding 
/o"8,Hn = 0.40±0.07 at a characteristic depth of z = 0.02. 
We note that the uncertainties on fas are marginalized 
over possible residual bulk flows. The results are also 
insensitive to details of, for example, corrections for ex- 
tinction in the SN hoso Further details about the un- 
certainties and light curve fitters used for the Al SNe 
can be found in THF. 

The second low- z MD result is from Davis et al. (2011 



hereafter DNM), who analyzed 2830 Tully-Fisher pecu- 
liar velocities at z < 0.033 and compared these to the 
predictions from the galaxy density field derived from 



the Two Mic ron All-Sky Survey Redshift Survey ( Huchra 
et al.||2012|). The analysis differs from THF: rather 



2012] . The analysis differs from THF 
than trie simple point-by-point comparison of predicted 
and observed peculiar velocities, DNM applied a spher- 
ical harmonic decomposition to both fields and found 
fas = 0.31 ± 0.05. We note that a straight average of 
the THF and DNM results yields fa s = 0.36 ± 0.04 in 
excellent agreement with the Composite vs. PSCZ result 
found above. 

We stress the that these two distance measurement re- 
sults are completely independent: they use different pe- 
culiar velocity samples, different density fields and dif- 
ferent reconstruction and comparison methods. The two 
measurements are consistent with each other: the differ- 
ence in fa s is 0.09 ± 0.08. The RSD and MD measure- 
ments are shown in Figure [l] 

3. DETERMINATION OF Q mfi AND <7g, 

In the ACDM model, 7 = 0.55. The model allows us 
to predict the value of f{z) = Q m (z)^ and as{z) at any 
redshift, assuming their values at redshift z = (denoted 
by subscript "0"). 



additi onal correction from the non-linear erg using Juszkicwicz et al. 
( 2010) and assuming f2 m .o = 0.27. 

3 Changing to Ry = 3.1 instea d of the default Ry = 1.7 changes 
fas by only 0.05. iKessler et all l |2009| l fit R v = 2.18 ± 0.14 sta t ± 
0.48 S yst, so the dirterence between Hy = 1.7 and Ry = 3.1 is a 
factor 2.8 times the systematic error. Hence the \-a Ry systematic 
effect on /<rg (~ 0.02) is much smaller than the random errors 
(0.07) 
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Figure 1. Growth parameter /<Xg(z) as a function of z. The 
data and errorbars are labelled as in Table 1. The ACDM model 
with WMAP7+BAO+i?o parameters f2 m ,o = 0.275, <r 8 ,o = 0.816 
is shown by the solid magenta curve. Note that the high-redshift 
RSD points assume the ACDM redshift-distance relation to correct 
for the AP effect and hence the appropriate value of fcrg(z). The 
dotted curves are normalized to have the correct eg at zqmb but 
have different possible 7 values starting at 0.50 at the top and 
increasing in steps of 0.05 towards the bottom. The dashed and 
dash-dot curves show predictions of flat and open DGP models, 
respectively, normalized to the amplitude of the fluctuations in the 
CMB (see text for details). 

For the fits, we use a simple x squared statistic with 
the following form: 



2 _ y^ [f(T S (meas)j - /o- 8 (pred) t 



(3) 



where the first term is the measured value and the sec- 
ond term is the model, and of is the uncertainty for each 
measured value. For RSD measurements at z > 0.1, /og 
is d egenerate with the redshift- distance relation, due to 
the |Alcock fc Paczynski] ( |1979[ hereafter AP) effect. If 
we assume a hat Universe, then the redshift-distance re- 
lation, and hence the AP distortion are all fixed by our 
choice of H m0 . With the AP distortion fixed, the ap- 
propriate fag values can then be calculated, but only for 
those surveys that have published the covariance matrix 
between the AP effect and /o 8 : WZ and BOSS. The pe- 
culiar velocity measurements (THF and DNM) and the 
low- z 6dF RSD measurement are negligibly affected by 
the AP distortion. 

Figure [2] shows the results of the fit in the f2 mj o ver- 
sus 0-3,0 parameter plane with 7 = 0.55. While the 
WZ+BOSS-only fit is quite degenerate in the fl m ,o — og.o 
plane, these degeneracies are broken when low- 2: data are 
included in the fit. Assuming a fixed 7 = 0.55, the best 
fit is Q mfi = 0.259 ± 0.045 and cr 8>0 = 0.748 ± 0.035 but 
the WMAP7+BAO+ifo value is consistent with these 
results. Some ACDM variants, such as those with a non- 
standard effective number of relativistic species, are dis- 
favored. The figure shows one such example (N e g = 3.8), 
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Constraints only from fa s 
measured from peculiar velocities 
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Figure 2. n m ,o versus <7g,o- Dashed contours show the 68, 95 
and 99 per cent confidence intervals using only WZ and BOSS 
redshift-space distortion measurements. The solid contours show 
the same constraints for the above data plus the peculiar veloc- 
ity measurements (THF and DNM) and the low-z 6dF RSD mea- 
surement (which is unaffected by the AP effect). The red circle 
shows the WMAP7+BAO+i?o best fit parameters for compari- 
son, the blue circle shows our best fit. The green circle shows the 
WMAP+SPT+ BAO+iJn best fit for a model with 3.8 effective rel- 
ativistic species ( jKeisler et al.|2011) , which is disfavored at > 99 
% confidence level 



from |Keisler et al. ( 2011 ). We stress that these results de- 
pend only - olT~gTowtrrmeasurements and hence are inde- 
pendent of other determinations (CMB, SNe, BAO etc.) 

4. CONSTRAINTS ON 7 

As discussed above, it is also interesting to measure 
the growth rate index 7. However, once 7 is included 
as a third parameter, the fits become very degenerate. 
These degeneracies can be broken using CMB data, for 
example from WMAP7+BAO+i^o- Let us assume that 
the Friedman equation (and hence the expansion his- 
tory) is the same as for the ACDM model, but treat 
7 as a phenomenological growth parameter which is al- 
lowed to differ from 0.55. At high redshifts (z ~ 1000), 
O m = 1 to high accuracy and so the growth of perturba- 
tions at early epochs is independent of 7. Therefore, we 
can use WMAP7+BAO+iJo parameters/M f2 mj o and o- 8 ,o 
to fix the amplitude of fluctuations at high redshift, as 
well as fixing the Friedman equation and expansion his- 
tory. Note that the quoted WMAP7 0-3,0 is extrapolated 
to z — assuming ACDM. To allow for other values of 7 
we use ACDM to calculate o- 8 (zcmb), and then extrap- 
olate the zcmb predictions forward at later times, u sing 
different values of 7 (following Samu shia et al.|2012 Sec- 
tion 4.5). We then use the WMAP7+BAO+±/ Monte 
Carlo Markov chains to marginalize over J7 m> o and o^o- 
The resulting fits are shown in Fig. [I] 

4 While the primary fluctuations in the CMB are independent of 
7, the secondary Integrated Sachs-Wolfe (ISW) effect is not. The 
ISW effect affects the CMB anisotropy spectrum on very large 
angular scales. We neglect this effect here. 
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Table 2 

Measurements of 7 from combinations of the 
data 



Sample 


7 


fCMB 


Ctot 


WZ 


666+ - 077 
u.ooo_Q.n 73 


0.053 


0.092 


LRG 


625+ - 072 

U.DZO_g Q77 


0.046 


0.088 


All RSD 


607+ - 038 
U ' DU ' -0.040 


0.046 


0.060 


THF+DNM 


U.DOd_ 064 


0.035 


0.077 


All 


619+° 033 
u.o±»_ 035 


0.042 


0.054 



1.0 



0.8 



0.6 



3NM 



^ 



f- 



0.4 



0.2 



LRG1 



l"HR 



6dF 



WZ3 



BOSS 



WZ2 



WZ1 



LRG2 



VVDS III WZ4 



2dF 



- - ACDM 
V77A best fit ±lcr 



).o 



0.2 



0.0 



0.8 



1.0 



Figure 3. Values of 7 derived for each survey individually, assum- 
ing the WMAP7+BAO+-ffo central parameters and extrapolating 
the WMAP7 crs^o to zcmb- Error bars on individual measurements 
reflect only the uncertainties in the fa% measurements and not the 
(correlated) uncertainties in the CMB-derived parameters. Sym- 
bols are as in Figure ffl The horizontal dashed line indicates the 
ACDM expectation. The best fit to the data and ±1<t range, in- 
cluding uncertainties in CMB-derived parameters is shown by the 
hatched box. 

Table 2 gives the derived 7 measurements for differ- 
ent combinations of the fag, measurements. Also listed 
are the uncertainties in 7 arising from fas, from the 
CMB-determinations of il m .o and <78,o, and the total er- 
ror. Note that the errors from WMAP+BAO+i^o are 
not independent (between one fa% measurement and an- 
other), are weakly dependent on z, and dominate the 
total error budget when all data are combined. Fig- 
ure [3] shows the derived 7 for all fas measurements as- 
suming the WMAP7+BAO+7?o parameters, and fixing 
C8(zcmb)- The results from the MD surveys are consis- 
tent with those from all of the RSD measurements. The 
hatched region shows the best fit to all datqj combined: 
0.619 ± 0.054, consistent with ACDM. 

Although for purposes of illustration we have focussed 
on a constant scale-independent 7, peculiar velocity data 
also allow one to test more complicated modified gravity 

5 If w e replace LRG1, LRG 2 and BOSS with the "free-growth" 
fits fromlTojeiro et a.1.1 (|2012|) , we obtain a 7 only 0.013 lower. 



scenarios (see e.g. the review by Clifton et al.|2012 |. One 
example is the DGP braneworld model. Fig. 1 shows the 
predicted values of fas for a self-accelerating fiat DGP 
model, with parameters chos en to match the expa nsion 
history and CMB constraints (Lombriser et al. 2009), and 
with a perturbation amplitude at early times chosen to 
match the CMB. Although the expansion history of the 
DGP model is similar to ACDM, the Friedmann equa- 
tion, and hence fl m (z), differs. Furthermore, its 7 = 0.68 
also differs considerably from the ACDM value. These 
effects lead to lower growth and much lower values of 
fas at z < 1. This model, already disfavor ed at the 5a 
level from other data ( Lombriser et al.|2009 |, is excluded 
at the ~ 8a level using only the fas measurements dis- 
cussed here. An open DGP model, which is a better fit 
to the redshift-distance relation, is excluded at the ~ 10er 
level by growth measurements alone. This suggests that 
alternative modified gravity models will be strongly con- 
strained by the requirement of simultaneously matching 
the ACDM expansion history and the ACDM growth his- 
tory. 

5. DISCUSSION 

The prospects for improving the measurements of 7 are 
excellent. At present, the error budget is dominated by 
the WMAP7+BAO+iJo estimates of the parameters at 
high redshift. Planck will reduce the CMB uncertain- 
ties so that these become subdominant. 

Peculiar velocity measurements will continue to im- 
prove, leading to a reduction in the observational errors 
in fas- The BOSS measurement will improve with fur- 
ther data releases. RSD measurements are also being 
made at higher redshifts (|Bielby et al. 2012). However 
the statistical power all of the RSD measurements com- 
bined (which are based on 800,000 redshifts) is similar 
to that of MD (based only on ~ 3000 peculiar velocity 
measurements), the subsamples having uncertainties of 
0.060 and 0.077 respectively in 7 (Table 2). It is there- 
fore clearly of great interest to improve the statistics of 
MDs. At low redshift, supernovae will continue to accu- 
mulate. We can also look forward to Fundamental Plane 



peculiar velocities from 6dFGS (Springob et al. 2012) 



and later an order of magnitude increase in the number 
of Tully-Fisher peculiar velocities from WALLAB\£j Fi- 
nally, the kinetic Sunyaev-Zel'dovich effect will be used 
to measu re the velocity fiel d of clusters at intermediate 



redshift (Hand et al. 



2012 ). The redshift survey data 
used to construct the density field and hen ce the pre- 
dicted peculia r velocities is also improving (Lavaux & 
Hudson|2011 1. It remains to better understand system- 



atics, by comparing measurements using the same sets of 
peculiar velocity and density data, but different methods. 

6. CONCLUSION 

We have shown that by combining measurements of 
fas(z) at different redshifts, and in particular by includ- 
ing results at z ~ from MD surveys, we can break the 
degeneracy between fl mt o and a$,o and obtain Sl m o = 
0.259 ± 0.045 and ct 8 , = 0.748 ± 0.035, consistent with 
independent determinations from WMAP7+BAO+-ffo- 

We can also constrain the growth index 7 by comparing 
measurements of fag(z) at low z, after fixing their values 

6 http://www.atnf.csiro.au/research/WALLABY/ 
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at zcmb ■ The strongest leverage on 7 arises from peculiar 
velocity measurements at the lowest redshifts. By includ- 
ing these measurements, we derive 7 = 0.619±0.054, con- 
sistent with ACDM. The Planck mission plus upcoming 
peculiar velocity and redshift surveys will tighten these 
constraints further. 
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